o
(N
Long range order in nonequilibrium interacting quantum spin chains
Tomaz Prosen^^^ and Marko Znidaric^
^ Department of Physics, Faculty of Mathematics and Physics, University of Ljubljana, Ljubljana, Slovenia
Department of Physics and Astronomy, University of Potsdam, Potsdam, Germany
(Dated: May 25, 2010)
We conjecture that nonequilibrium boundary conditions generically trigger long range order in
nonequilibrium steady states of locally interacting quantum chains. Our result is based on large
scale density matrix renormalization group simulations of several models of quantum spin 1/2 chains
which are driven far from equilibrium by coupling to a pair of unequal Lindblad reservoirs attached
locally to the ends of the chain. In particular, we find a phase transition from exponentially decaying
to long range spinspin correlations in integrable Heisenberg XXZ chain by changing the anisotropy
parameter. Long range order also typically emerges after breaking the integrability of the model.
PACS numbers: 05.70.Ln, 05.30.Rt, 03.65.Yz, 75.10.Pq
(N
I
o
>
o
^,
o
o
X
Nonequilibrium physics of low dimensional locally in
teracting quantum manybody systems is becoming an
increasingly hot topic, in particular due to enormous re
cent progress in experimental techniques with ultracold
atomic gases and optical lattices [l , and arguable connec
tion to important open challenges in condensed matter
physics such as hightemperature superconductivity and
solidstate quantum computation. While thermalization
in isolated manybody systems seems to be well under
stood 13 , and near equilibrium physics is essentially given
by linear response formalism, far from equilibrium even
some very basic concepts remain unclear.
One of the key questions that we address in this Let
ter is: Driving a large onedimensional locally interacting
quantum system out of equilibrium by coupling to two
unequal reservoirs at the system's ends, will the non
equilibrium steady state (NESS) in the bulk be well de
scribed by the notion of local equilibrium? If yes, then the
properties of NESS are essentially given by local Gibbs
state with a well defined local temperature (and perhaps
local chemical potential) field(s). If not, then the local
temperature has no meaning, and the physics of far from
equilibrium system is qualitatively different than in equi
librium. While MerminWagner's theorem [3] prohibits
certain types of order at finitetemperature equilibrium,
(arbitrarily weak) nonequilibrium driving at the bound
ary of the system breaks not only the translational invari
ance but typically also the symmetries of the Hamiltonian
and opens the possibility of long range order (LRO). As
has been shown for an open XY spin 1/2 chain, LRO
may set in abruptly by varying the system's parameters
through a non equilibrium quantum phase transition [j.
Nevertheless, since the XY chain is equivalent to a
quasifree fermionic model, one might argue that LRO
might be a nongeneric effect which would typically not
survive switching on the interaction, or even more so,
breaking the integrability of the model. In this Letter we
show that this is not the case by performing large scale
densitymatrixrenormalizationgroup (DMRG) simula
tions of NESS for integrable and nonintegrable versions
of anisotropic locally interacting Heisenberg (XXZ) spin
1/2 chain. By changing the anisotropy parameter we find
a clear evidence of nonequilibrium quantum phase tran
sition from exponentially decaying spinspin correlations
to LRO where spinspin correlations saturate and do not
depend on the system size. Perhaps even more strik
ingly, we find that LRO is typical for nonintegrable de
formations of the model. Note that the physics of NESS
can be very different from the much studied quench dy
namics where one can observe at most a transient non
equilibrium phenomena.
The model we study is an anisotropic nearestneighbor
interacting spin 1/2 Heisenberg XXZ chain
^ = EK^i+JJ+i+^K+i)
(1)
The Hamiltonian ([T]) is an important model for several
reasons. It describes various spinchain materials [5] like
cuprates and as a particularly promising experimental
technique, it can be realized with cold atoms in an opti
cal lattice [6. From the theoretical side, ([T]) is perhaps
the simplest model of strongly interacting electrons (to
which it can be mapped by Jordan Wigner transforma
tion). Understanding the Heisenberg model would there
fore provide a guidance to physics of more realistic mod
els of interacting electrons like e.g. the Hubbard model.
Despite integrability of the Heisenberg model J7 , cal
culating quantities beyond, for instance, the ground state
energy or lowenergy excitations is exceedingly difficult.
Spinspin correlation functions are particularly impor
tant quantities. Regarding equilibrium correlation func
tions in the Heisenberg model the following is known
7]: (i) at infinite temperature there are no correlations,
(ii) at finite temperature correlations asymptotically de
cay exponentially with distance, (iii) at zero tempera
ture (ground state) correlations decay asymptotically as a
powerlaw of the distance in the gapless phase for  A < 1,
while they decay exponentially in the gapped phase of
A > 1. For nonequilibrium stationary states almost
nothing is known about the correlation functions, some
preUminary hints can be found in [8. The aim of this
work is to focus on long range correlations in NESS.
While the properties in closetoequilibrium situation
can be studied using a linearresponse formalism, far from
equilibrium one is forced to consider genuine nonequilib
rium setting where the central system is coupled to reser
voirs at different potentials, for instance, at the chain
ends. We are going to take into account the Markovian
reservoirs in an effective way by using the Lindblad equa
tion. Time evolution of systems's density matrix is there
fore governed by the Lindblad master equation [9
^,pAp.h]
•bath
{p)
(2)
The dissipator £!°^^^ — £}^ ^ £,^ is a sum of two parts,
one acting on the left end and the other acting on
the right end of the chain. Each one is written in
terms of a set of Lindblad operators L^ as C^{p) =
T.,{[LlP;Ll']
[LlpLf:
X
L, R. Most of the
time we shall use the onespin bath for which Lindblad
operators act only on a single spin. Two Lindblad opera
tors at each end are i^2 = \/l ^±^1 ^ ^?
7 i±c
n
where a/, = h{(T^ ± ia^). In numeri
l±tanh£f^
cal simulations we always consider strong coupling 7=1.
These operators are such the the stationary state of the
(say, left) bath spin only, C^{pi) = 0, is pi oc exp (fJhcrl)
with polarization (cr^) — tanh/iL. One can argue that
such a situation formally corresponds to an infinite tem
perature with fixed bath magnetizations. Later we shall
also consider simple twospin baths where two border
spins are coupled to a bath in order to simulate the effect
of finite temperature. The choice of the bath is qualita
tively inessential for the results shown.
To find NESS (fixed point of the Liouvillean) we sim
ulate the Lindblad equation ^ using timedependent
DMRG (tDMRG) [iBl with matrix product operator
(MPO) ansatz. This enables us to study systems much
larger than it would be possible with bruteforce diago
nalization. After sufficiently long time, the density ma
trix p{t) converges to a timeindependent NESS. Once
NESS density matrix pness is obtained in terms of MPO,
various expectations can be efficiently evaluated. [l8
We now turn to the study of the main object of the
present work, the 2point spinspin correlation function
C{i,j)^ trpNESSojOj  (tr/9NESsCTD(trpNESS0j) (3)
C{i,j) is usually found to depend on the relative dis
tance r = i ~ j ranch stronger than on the perpendic
ular coordinate, thus we often show a onedimensional
crosssection of C{i,j) along the skewdiagonal, that is
C(r) := C(fc, k + r), where k = [s±ij  [§J . Sometimes,
when larger fluctuations among the values of C{i,j) are
10'
10
:10
O
10
10
10
10
\ (a)
exp(r/;,)
1
\
z
2
\

■A
^*o^^
n=6 "'Sl!>^
o" """'^■■■■. ^^\
4
^^^^
n=10 X
n=60\^ ;:::::„;■■■..
n=20 ,
(
] 2 4
6 8 10 12 14
16
" I
r
:: (b)
(A^if(^r/nf
— ;
1
1
2
ii.  VvVVv 
ttt:::^;^^
n=lD..,
•20:: : ::i^^
3
■ "■^r'"'""'^v"
, ■ :
4
\v
010
10
10
0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1
r/n
FIG. 1: Spinspin correlation function C{r) in open XXZ
model at A = 0.5 (a) and A — 1.5 (b), both for fiL.n =
±0.22. Symbols are from exact calculation, dashed curves
from tDMRG, while full curves suggest n — >■ 00 asymptotics.
noticed, we average the absolute value of the correla
tion over all pairs of points with a constant distance
(C(r)) := {\C{k, k \ r)\)k, and in order to characterize
LRO, as order parameter we compute the residual cor
relator as an aditional average over r between n/4 and
3n/4, (q):=(C(fc,fcHr))^/^^'■=^^"/^
It turns out that the spinspin correlation function in
NESS depends strongly on the anisotropy A. For exam
ple, for A = 0.5 the correlation function (its cross sec
tion C{r) is shown in Fig. [1^) shows exponential decay,
with asymptotic scaling C{r) ^ e"''/^, with ^ = 2.0. We
can see that the correlation decay rate ^ is independent
of the chain size. Therefore, correlations in NESS are
shortranged. However, after increasing the anisotropy
the situation becomes dramatically different, for exam
ple, for A = 1.5 the correlator C{r) clearly exhibits LRO
as shown in Fig. [T]d. With increasing the chain size n
the correlation function converges to the limiting form
for which we can observe several things: (i) spatial de
pendence is a function of the scaled variable r/n only,
for large n, meaning that there is LRO in NESS, (ii)
the size of the correlations scales with the driving as
^ (Ap)^ = (^L — Mr)^ and is independent of n. Further
more, correlations are positive (the point where they be
come negative moves towards r/n — ;■ 1 as n is increased).
10 15 20 25
FIG. 2: Residual correlator (C) in open XXZ model against
the scaling variable (A — Ac)n with Ac = 0.91, for n = 20
(circles), 30 (diamonds), 40 (squares). The inset shows un
sealed data. Here /^l.r = ±0.22 whereas practically the same
scaling of {\C\)/{Afi)'^ was obtained also for /il.r = ±0.05.
The limiting form of the correlation function is perhaps
close to C{i = nx,j = ny) ~ (A/i)^a;(l — y). Very similar
form of the 2point correlation function, however with a
size dependent prefactor (l/n), has been found in cer
tain noninteracting, dissipative (noisy) or even classical
models: in a ballistic exactly solvable XY model 4,], nu
merically observed in the diffusive Heisenberg model with
dephasing 11] and analytically explained in a diffusive
XX model with dephasing [12 , but also in some classical
lattice gasses like simple exclusion processes 13 or its
quantum analogs 141]. The fact that long range corre
lator here does not decay with n perhaps hints at truly
quantum coherent character of these correlations.
From Fig IT] we can conclude that somewhere between
A = 0.5 and A ~ 1.5 the range of correlations must
change from short to long range, i.e., there is possibility
of a nonequilibrium phase transition. To find if this is
the case and where is the transition point we study the
residual correlator (C) against parameter A for several
system sizes. Remarkably, numerics (shown in Fig. [2])
indicates the scaling (C) ^ /((A — Ac)n''), where / is
some universal scaling function and t/ « 1.0, suggesting
that the transition becomes noncontinuous in the ther
modynamic limit. Results seem to be accurate enough to
support the claim that the observed critical anisotropy
Ac « 0.91 is distinct from the zero temperature critical
anisotropy A* = 1 where the ground state gap opens.
One may argue that LRO observed above might be due
to integrability of the XXZ model. In order to check the
effect of (non)integrability we add a staggered magnetic
field to the Hamiltonian Q, namely H' = H+X]?=i ^j"'?
with h2j = h, h2j+i = 0, for h = —0.5, which is known to
bring the model in the regime of quantum chaos [8 . Us
ing the same driving as before, we find that NESS again
exhibits long range correlations. In order to somewhat
10"
10^
1 102
4
■
w

— <:"""■
n=20\
•,30
•^■:40
,,.50
10
5 10 15 20 25 30 35 40 45 50
r
FIG. 3: Smoothed spinspin correlation function (C(r)) in
the nonintegrable deformed XXZ model with A — 0.5 and
staggered transvere magnetic field h = —0.5, /iL,R = ±0.1.
isolate the effect of integrability breaking alone, we take
A = 0.5, which in the integrable limit exhibited short
range correlations. In order to smoothen the oscilla
tions in the correlation function due to staggered field we
plot in Fig. [3] the average correlator (C(r)) and clearly
observe LRO which is again consistent with a constant
residual correlator (C) — 0{n^).
All results presented so far have been for symmetric
driving, fj,^ = — /zl, which results in NESS having aver
age magnetization zero and the average energy density
also zero. This can be interpreted as NESS at an infi
nite average temperature, although care must be taken
when discussing the temperature in an integrable open
system because thermalization is absent for unadjusted
reservoirs due to lack of ergodicity [l5 . Thus the struc
ture of the ground state very likely plays no role in the
above NESS correlations. An interesting question is how
do the correlations in NESS change as we lower the ef
fective temperature? Also, correlations could depend on
the average magnetization, as does for instance the spin
transport in a linear response regime. In the zero mag
netization sector there appears to be a transition from
ballistic spin transport for A < 1 to a diffusive trans
port in the gapped regime A > 1. On the other hand,
for nonzero magnetization spin transport is always bal
listic as exemplified by a nonzero Drude weight 16. To
study correlation function in NESS at a nonzero energy
density (i.e., at finite "temperature") we used C^'^^^ cou
pling to two spins at the boundary or a onespin bath
with nonsymmetric driving. Similar results are obtained
in both cases and we are going to show only the ones
for a twospin bath as the relaxation to NESS there is
faster. Lindblad operators on two boundary spins are
chosen in such way that the invariant state of these two
spins would be a grandcanonical state at some target
temperature and inagnetization in the absence of the
bulk Hamiltonian [ll, l7. Effective temperature field
of NESS is then prescribed by equating the measured
local magnetization and energy density in the bulk to
FIG. 4: Finite temperature correlation function C{r) for open
XXZ model at A = 1.5. Frame (a) is for T^eas « 33 while (b)
is for Tmoas ~ 5.8 (both measured in the middle of the chain).
At short distances (r < 10 in (b)) we have an exponential
decay of C{r) with the sizeindependent decay rate. For large
r, we find the residual correlator scaling as (C) ~ (Ap) /n"
(insets), with a ~ 0.4 in (a) and a « 2 in (b).
the grandcanonical one at some "measured" tempera
ture Tnicas and potential /imoas through canonical den
sity matrix p ~ exp [—{H — fj,S^)/T], see [15 for details.
The results for correlation function at two different mea
sured temperatures in the center of the chain are shown
in Fig.Ufor the gapped XXZ model A = 1.5 at zero stag
gered field h = 0. These two results for finite "tempera
ture" show that the decay of the correlation function has
two regimes. For small distances r, we have a shortrange
exponential decay, with the decay rate depending only
on the temperature and not the chain length. This can
be interpreted as a finitetemperature equilibrium part of
the correlation function. On the other hand, at larger
distances longrange correlations still persist, but scaling
algebraically with the system size as (C) ^ (A/i)^/n",
with some nonzero power a. This second regime can be
interpreted as being genuinely nonequilibrium. Noting
that at an infinite temperature. Fig. [Hd, we have a =
as discussed before, while at lower temperatures, Fig. 4l
a increases. We conjecture that a goes from at an
infinite temperature to a — >■ oo at zero temperature. In
creasing exponent a with decreasing temperature might
be in line with the exponential decay of correlations for a
gapped ground state in equilibrium. Note however, that
in the equalbath (equilibriumlike) situation, A/i — >■ 0,
all correlations in NESS decay exponentially.
In conclusion, we have presented strong numerical sup
port for the existence of long range correlations in bound
ary driven nonequilibrium steady states of locally inter
acting quantum spin chains. We have shown the exis
tence of quantum phase transition in an open XXZ spin
1/2 chain with spontaneous emergence of long range cor
relations with increasing the anisotropy, and/or breaking
the integrability of the model. These long range cor
relations are particularly strong, as they  evaluated at
the fixed scaled distance r/n  do not depend in the sys
tem size n, and are thus fundamentally different from the
'hydrodynamic' long range correlations observed for in
stance in certain classical far from equilibrium stochastic
processes [13j . This phenomenon could have far reach
ing consequences for far from equilibrium physics in low
dimensional quantum manybody systems, such as com
plete absence of local equilibrium and nonexistence of
local thermodynamic potentials. The work is supported
by grants Pl0044, Jl2208 of Slovenian Research agency.
[1] I. Bloch et ai, Rev. Mod. Phys. 80, 885 (2008).
[2] M. Rigol et ai, Nature 452, 854 (2008).
[3] N. D. Mermin and H. Wagner, Phys. Rev. Lett. 17, 1133
(1966); P. C. Hohenberg, Phys. Rev. 158, 383 (1967).
[4] T. Prosen and I. Pizorn, Phys. Rev. Lett. 101, 105701
(2008); T. Prosen and B. Zunkovic, New J. Phys. 12,
025016 (2010).
[5] C. Hess, Eur. Phys. J. Spec. Top. 151, 73 (2007); A.V.
Sologubenko et ai, J. Low Temp. Phys. 147, 387 (2007).
[6] S. Trotzky et ai, Science 319, 295 (2008).
[7] A. Luther and I. Peschel, Phys. Rev. B 12, 3908 (1975);
V. E. Korepin, N. M. Bogoliubov and A. G. Izergin,
Quantum inverse scattering and correlation functions,
Cambridge University Press, 1993.
[8] G. Benenti et al, Phys. Rev. B 80, 035110 (2009).
[9] R. Alicki and K. Lendi, Quantum dynamical semi
groups and applications, Springer, 2007; H.P. Breuer and
F. Petruccione, The theory of open quantum systems, Ox
ford University Press, 2002.
[10] A. J. Daley et ai, J. Stat. Mech. P04005 (2004).
[11] M. Znidaric, New J. Phys. 12, 043001 (2010).
[12] M. Znidaric, J. Stat. Mech. L05002 (2010).
[13] B. Derrida, J. Stat. Mech. P07023 (2007).
[14] K. Temme et ai, arXiv:0912.0858
[15] M. Znidaric et al, arXiv:0910.1075
[16] X. Zotos et al., Phys. Rev. B 55, 11029 (1997).
[17] T. Prosen and M. Znidaric, J. Stat. Mech. P02035 (2009).
[18] In order to make sure that the truncation errors
in tDMRG algorithm are under control, we carefully
checked our tDMRG results against the result of an ex
act calculation for small systems up to n < 10. In each
tDMRG run we also ensured that the results converged
with respect to increasing MPO dimension D (in some
cases, we went up to D = 150) and time of simulation.